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Abstract 

We use a configuration space chiral model in order to evaluate nucleon and delta a— terms. 
Analytic expressions are consistent with chiral counting rules and give rise to expected non-analytic 
terms in the chiral limit. We obtain the results aj\f = 46 MeV and a a = 32 MeV, which are very 
close to values extracted from experiment and produced by other groups. 
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I. INTRODUCTION 



The delta (A) plays a very important role in low-energy pion-nucleon (ttN) scattering 
and correlated processes, such as the nucleon-nucleon interaction. Its contribution as an 
intermediate state, in many instances, supersedes that of the the nucleon. This happens for 
two main reasons. The first one is that the 7riVA coupling constant is rather large, whereas 
the other is related to chiral symmetry. 

At low energies, pion-hadron interactions are well described by effective theories, in which 
an approximate SU(2) x SU(2) symmetry, broken by the pion mass (p), accounts for the 
smallness of the u and d quark masses. In this framework, elastic pion-baryon scattering 
is dominated by diagrams involving both contact terms and propagating states. In order 
to comply with threshold chiral theorems, the latter are typically given by polynomials in 
small quantities, such as the pion mass or three-momenta, divided by energy denominators. 
When the delta is present, the scale of some denominators is given by the quantity u>a = 
(M 2 — m 2 — /i 2 )/2m, where M and m are respectively the delta and nucleon masses. As the 
difference A = M — m is small, one has cg>a ~ A. Delta contributions are given by ratios of 
small quantities and may turn out to be large. In such cases, numerical values adopted for 
A do influence predictions produced by effective theories, especially those that rely on the 
small scale expansion[l] or the heavy baryon approximation [2]. 

In chiral perturbation theory, there is a clear conceptual distinction between the bare 
baryon masses, present in the lagrangian, and their observed values, which include loop 
corrections. The former should, in principle, be preferred as inputs in the evaluation of the- 
oretical amplitudes. Nevertheless, as there is little knowledge available concerning the bare 
delta mass, one tends to use physical values in calculations. In most cases, it is reasonable 
to expect that this would have little numerical importance. On the other hand, in the case 
of the parameter A, which is a small quantity, the influence of loops may become relatively 
large. 

Recently Bernard, Hemmert and Meissner[3], have stressed that the value of Ao, the 
delta-nucleon mass splitting in the chiral limit, is an important constraint to lattice data 
extrapolation. The purpose of the present work is to estimate the delta er-term, which 
controls the change induced in A when one goes from bare to physical masses. This er- 
term was studied in the framework of a quark model by Lyubovitskij, Gutsche, Faessler and 
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Drukarev[4] and the reader is referred to their paper for a clear formulation of the problem 
and earlier works. 

According to the Feynman-Hellmann theorem [5] the mass mg of a baryon B is related 
to its sigma-term <7b by ob = fi 2 dms/df! 2 . Therefore the sigma-term provides a measure 
of the shift in the baryon mass due to chiral symmetry breaking. Whenever it is possible to 
evaluate <jb as a function of fi, the bare mass mg, can be extracted from the relation 



As the leading term in <jb is proportional to fi 2 , it enters directly the mass shift and 
the difference m# — <7b already provides a crude estimate for the bare value. In the case 
of the nucleon, one has (7^=45 MeV[6], which amounts to 5% of its physical mass. In 
chiral perturbation theory, the leading contribution to cannot be predicted theoretically. 
Formally, it is associated with the constant C\ of the second order lagrangian[7, 8], which can 
be extracted from empirical subthreshold information. The situation of the delta is much 
worse, for 7rA scattering data are not available. One is then forced to resort to models in 
order to calculate the delta a-term, which is associated with the parameter a\ defined in 



In this work we estimate a a using a model which proved to be successful in the case of 
the nucleon. Our paper is organized as follows. In section II we review our calculational 
procedure in the case of the nucleon and present results for the delta in section III, leaving 
technical details to the appendices. The main expressions for both the nucleon and delta 
a-terms in configuration space are given in appendix B, written in terms of the loop integrals 
defined in appendix A. The consistency of our results with standard chiral counting rules is 
discussed in appendix C whereas their behaviour in the chiral limit is given in appendix D. 
A summary is provided in section IV. 

II. MODEL FOR THE SIGMA TERM 

In order to evaluate a a, we follow a procedure used previously in the study of <Tjv(t), the 
nucleon scalar form factor [9], which is briefly reviewed here. The leading contributions to this 
function is 0(q 2 ) whereas the triangle diagram, involving only known masses and coupling 
constants, gives rise to corrections which begin at 0(q 3 ) and are completely determined. At 




(1) 



ref.[l]. 
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0(q 4 ), on the other hand, interactions incorporate the low energy constants C\,C2 and C3. 
Data on irN subthreshold coefficients indicate that C2 and C3 are larger than c\ and that their 
values are approximately saturated by A intermendiate states[8]. Thus, up to C(q' 4 ), the 
function a^ii) can be well represented by the leading tree contribution associated with c±, 
supplemented by the two triangle diagrams shown in fig. 1, involving N and A intermediate 
states. In the sequence we will make use of the fact that, in configuration space, contact and 
loop contributions split apart, since the Fourier transform acts as a filter [10]. As a result, 
the theoretically undetermined leading tree term yields a zero-range 5-function, whereas the 
triangle diagrams give rise to spatially distributed structures, fully determined by known 
parameters. 




FIG. 1: Contact term and triangle diagrams contributing to the cr-term. 



The nucleon scalar form factor in momentum space is defined by 

(N(p')\-jr sb \N(p)) = a N (t) u(jp') u(p) , (2) 

where C s b is the symmetry breaking term in the lagrangian and t = (p'—p) 2 . In terms of 
the quark degrees of freedom, one has C s b = —rh(uu+dd), with rh = (m u + ra<i)/2. The 
configuration space scalar form factor is denoted by &n and given by 

Mr) = J ^f^ r a N {t) (3) 

with q = (p'—p), in the Breit frame. The nucleon a-term, defined as oat = cjv(i = 0), is 
given by 

POO 

a N = 47r / dr r 2 a N {r) . (4) 
Jo 

The contributions from the diagrams of fig. 1 to crjv(r) read 

vn{t) = -4ci/i 2 5 3 (r) +a NN (r) + 5"jv A (r) , (5) 
4 



where 07^ (r) and 5"jv A ( r ) are given by eqs.(B4) and (B5) of appendix B and displayed in 
fig. 2. These functions are based on unregularized loop integrals and diverge for small values 
of r. In momentum space, regularization is achieved by means of counterterms, which 
give rise to polynomials in t, designed to cancel the divergences of the loop integrals. In 
configuration space, this regularization procedure amounts to adding (^-functions and their 
derivatives to crjv(r). This gives rise to a regularized form factor which is very large both at 
r = and in a sizeable vicinity of that point. We argue, in the sequence, that this picture 
is not consistent with the definition of the form factor given by eq.(2). 




r(fm) 

FIG. 2: Spatial dependence of the nucleon scalar form factor (continuous line) and partial con- 
tributions due to N and A intermediate states [eq.(B4), dashed line and eq.(B5), dotted line, 
respectively]. 



Pions are Goldstone bosons, collective states derived from the qq condensate. The corre- 
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sponding degrees of freedom are appropriately accomodated into non-linear lagrangians and 
described by the field U = exp(zr - it 9), where 7r is the isospin direction and 6 is the chiral 
angle. This function can be expressed as U = cos 6 + it ■ tt sin# and the dimensional pion 
field is given by 

<p = f-n sin 6 7r . (6) 

Long ago, Skyrmefll], in a series of papers, considered the possibility of pion fields being 
either weak or strong. It is worth noting that the words weak and strong, as used here, 
are akin to the notion of weak and strong electromagnetic fields developed by Schwinger, 
and not at all related to the nature of the fundamental interactions. In the former case, 
changes in the qq condensate are small, one relies on the approximation </> ~ f n 6 tt and 
can employ perturbative techniques, as in chiral perturbation theory (ChPT). In the latter, 
disturbances of the QCD vacuum become important and the non-linear nature of pionic 
interactions manifests itself through the condition \4>\ < f n . The physical picture behind 
eq.(6) is that pions, as Goldstone bosons, destroy the qq condensate in order to exist. 

When strong fields are present, constraints also apply to the scalar form factor. The 
symmetry breaking lagrangian is written in terms of the dimensional pion field as 

C sb = \fl f Tr [U + 17+] = fl /i 2 cos 9 . (7) 

This structure shows that £ s t, is a bound function and definition (2) means that the same 
necessarily happens with the scalar form factor. The function a^ix) corresponds to a mass 
density induced in the vacuum by the presence of the nucleon, which manifests itself in the 
form of a pion cloud. Far away from the nucleon, eq.(7) yields the density of the condensate, 
which is negative and equal to — fl fi 2 . In the description of a nucleon, it is convenient to 
use a convention for the energy in which the density tends to zero at long distances and C s b 
is rewritten as 

£ sb = f^ 2 (cos6-l) . (8) 

In this new convention, the density vanishes when r — > oo and increases monotonically as 
one approaches the center of the nucleon as in fig. 2. At a critical radius R one has cos 6 = 1, 
the density becomes that of empty space and the condition 

a N (R) = f 2 ^ (9) 
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holds. Beyond this point, a further increase in <jn would correspond to cos 9 > 1. In order to 
prevent this behaviour, we assume that the condensate no longer exists in the region r < R, 
and that the energy density saturates at r = R. For this reason, in our previous evaluation 
of cr/v[9], we used the expression 

a N = -tiR 3 fy + 4vr / dr r 2 a N (r) (10) 
<J Jr 

instead of eq.(4). This procedure is the basis of our model. 

In the numerical determination of ctjv, we use the results of appendix B and consider two 

possibilities for the 7riVA coupling constant in the lagrangian (B2), corresponding to either 

the SU(A) prediction g^NA = SgA/^V^ = 1-33 or g^NA = 1.47, which yields r=120MeV for 

the A decay width. The corresponding results, given in table I, are quite close to the value 

extracted from experiment by Gasser, Leutwyler and Sainio[6], namely = 45MeV. 

TABLE I: Nucleon u-term as function of the ttN A coupling constant. 



9-kNA 


R (fm) 


a N (MeV) 


1.33 


0.57 


45.8 


1.47 


0.59 


49.4 



Consistency with chiral symmetry is an important issue in this problem. Therefore, we 
note that, although the chiral powers of the pion mass expected from triangle diagrams 
are not explicit in the expressions of appendix B, the use of covariant relations among 
integrals[12] allows results for partial contributions to the a-term to be recast in such a way 
that these powers become apparent, as shown in appendix C. In appendix D we show that 
the formal chiral expansion of eq.(10) gives rise to the expected non-analytic terms (log /i and 
/i 3 ) and agrees fully with that produced by standard chiral perturbation theory [8], provided 
the renormalization scale is identified with 1/R. 

III. DELTA cj-TERM 

The delta scalar form factor is defined as 

< A(p',s')\ - £ sb \A(p,s) >= -<(p') ^V A (t)+pVF T (t)] <(p) , (11) 
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where u s u is the A spinor[13] and o"a and Ft are respectively the scalar and tensor form 
factors. The minus sign on the r.h.s. is associated with the conventions used in the free 
A lagrangian as in ref.[l]. We assume that the scalar form factor is determined by a short 
range contact interaction and the two long range two-pion processes shown in fig. 1. 

In figs. 3 and 4 (zoom in) we display the profile functions for the partial contributions to 
a a given by eqs.(B9,B10) and it is interesting to note that the nucleon contribution oscillates 
in the outer region, in sharp contrast with fig. 2. This behaviour is due to the fact that the 
delta is unstable and makes a a to be smaller than a^. The structure of partial contributions 
for SU(4) coupling constants is given in table II, where core and cloud refer respectively to 
regions inside and outside the cutting radius R. 




0.2 0.4 0.6 0.8 1 1.2 1.4 

r(fm) 



FIG. 3: Spatial dependence of the delta scalar form factor (continuous line) and partial con- 
tributions due to N and A intermediate states [eq.(B9), dashed line and eq.(BlO), dotted line, 
respectively]. 
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FIG. 4: Expanded portion of fig. 3, with the same conventions. 



TABLE II: Partial contributions to a N [eqs.(B4,B5)] and a A [eqs.(B9,B10)]. 





core 


cloud N 


cloud A 


sum 


a N (MeV) 


16.7 


13.0 


16.1 


45.8 


a a (MeV) 


14.3 


-1.5 


19.3 


32.1 



Processes containing nucleon intermediate states give rise to an imaginary component a A 
for the delta cx-term, which can be related to the decay width by means of the Feynman- 
Hellmann theorem [5]: 



2 d{M-iV/2) 



aA~ia A = » ^ • (12) 



9 



Using [8] 



9 " NA &[{M+m) 2 -tf\ 



24ttM 2 /^ 

qA = 2M ^ MA+mA+ ^- 2m2M2 - 2 ^ 2M2 - 2 ^ m2 ' ( 13 ) 



one finds 



°^-MrA& + ^^r^ 2 -^} • (14) 

The values of the distance R for which cr&(R) / f^fi 2 = 1 and of the delta cr-term, calculated 
by means of eqs.(B9, BIO), are given in table III, for different choices of the coupling 
constants g n NA and g n AA- The S77(4) predictions for these constants are 1.33 and 0.75, 
whereas the value 1.47 for the former yields the empirical decay width. The value 0.67 for 
the latter was used in ref. [3]. Results for the real component of a a are sensitive to the 
coupling constant g^AA and fully consistent with that given in ref. [4], namely a a = (32 ± 3) 
MeV. On the other hand, our prediction is larger than that quoted in ref. [3]. The values for 
the imaginary component a A , obtained by means of eq.(B9), are identical with those given 
by eq.(14), as they should. 



TABLE III: Real and imaginary parts of the A <r-term as function of the coupling constants. 



9-kNA 


9nAA 


R (fm) 


Re a a (MeV) 


Im a a (MeV) 


1.33 


0.75 


0.54 


32.1 


-21.7 


1.33 


0.67 


0.52 


28.1 


-21.7 


1.47 


0.75 


0.55 


31.7 


-26.7 


1.47 


0.67 


0.53 


27.8 


-26.7 



IV. SUMMARY 

We have discussed a model aimed at determining a— terms, which consists in cutting off 
configuration space expressions at the point where the cosine of the chiral angle becomes 
larger than 1. The model has been used to calculate <jn and a a with success. In the 
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former case, a value very close to that extracted from experiment by Gasser, Leutwyler and 
Sainio[6], was obtained. In the case of the delta, the prediction 28 MeV < <ta < 32 MeV, 
depending on the coupling constants employed, is also very close to the result produced 
by another group [4]. The fact that the delta can decay gives rise to a pion cloud which 
includes an oscillation and is responsible for both the relation a a < &n an d the consistency 
of the imaginary part of a a with the decay width. Analytic expressions also comply with 
chiral counting rules and give rise to expected non-analytic terms in the chiral limit. These 
features suggest that our calculational procedure is sound and can be reliably applied to 
other systems. 



APPENDIX A: LOOP INTEGRALS 



In the triangle diagrams, p and p' are the initial and final baryon momenta, whereas k 
and k! are the momenta of the exchanged pions. We also employ the variables 

q=(p-p'), P = (p+p')/2, Q = (k + k')/2. (Al) 

In all diagrams, the external baryon, with mass m e , is assumed to be on shell and one 
has 

p 2 = p ' 2 = m 2 e , p ■ q = , (A2) 

u yfu = , u P u = muu , (A3) 

u»tiu v = 0, rt L Pu v = Mu»u v . (A4) 

The basic loop integrals needed in this work involve either two or three denominators. 
We use the definition 

d ' Q 1 (A5) 



(2ny [(Q+q/2) 2 -^} [(Q-q/2) 2 -^} 
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and the dimensionless expressions 

I 7T7T ~ 



: Tim 

(4tt) 2 ™ ' 



TfJtV 
7T7T 



X7T7T 



/i 2 (4vr) 2 

2\im e i 
m + P) 2_ ml] = (4^ 



[<T n^ 0) + • • •] 



_ n (ooo) 

2 X7T7T ' 



[( g +P) 2_ m 2 ] (47r) 



n( 10 °) j 

J -- L :r!7r7r 



[(g +P) 2_ m 2 ] (47r) 



[<T ni™ 0) + • • •] 



M^^!. fK 10 °) 4 



(A6) 
(A7) 
(A8) 
(A9) 
(A10) 
(All) 



J [(Q + P) 2 -m 2 ] (4vr) 2 

where the ellipses indicate terms that do not contribute to the scalar form factors. The 
usual Feynman techniques for loop integration allow one to write the regular parts of these 
integrals as 



7T7V 

fr(° ) = 

7T7T 

n (fc00) 

xnn 







1 fD^ 
da In I — — 

1 D ( D 
da - — - In 



2/i 2 



da a f db [— m e (l— a)/fi] k 
Jo 

i 



— - / da a I db \— m e (l— a)/u] k In ( — — 
A* Jo Jo V2/im f 



(A12) 
(A13) 
(A14) 
(A15) 



with 



D n7V = ji 2 — a{\ — a) q z 



a jj, + (1 — a) m x — a (1 — a) m e — a b(l—b)q 



(A16) 
(A17) 



The dimensionless configuration space functions S are defined as 

S = 



e -ik-x n _ 



(27T) 



(A18) 
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with x = fir and k = q/fi. Performing the Fourier transforms, we find 

SlT = — 2 Ki(2x) , (A19) 

TlX 

SlT = --^ [xK (2x) + K^x)] , (A20) 

- 1 \ ™ (^ „\ /..]k 



fi TlX Jq 



jj x J a 



2 > - S2£> = ~ — - l da [ - me(1 ~ a)//i] K^x) , (A21) 

H Tix Jq a 

[Fo(2|0|x) + zJ o (2|0|x)] , (A22) 

2 > - C° 0) = A — r rfa «h^e(l-a)//i] fc /^i(20x) , (A23) 

TlX 1 fj, J 

2 <O-C°. O) = -^^daai-m^l-a)/^ |0| [F 1 (2|0|x) + zJ 1 (2|0|x)] ,(A24) 

with 

2 = [ a/i 2 + (l-a)m^ -a(l-a)m^]/(/i 2 a 2 ) . (A25) 
APPENDIX B: SCALAR FORM FACTORS 

We give here the expressions for a^, due to triangle diagrams containing external states 
B and an intermediate states /. The following interaction lagrangians[l, 8, 13, 14] are used 

Cknn = |f {NwsT a N} , (Bl) 

^ J IT 

C„ NA = ^ {A M - (Z- 1/2) 7m7i/ ] M a A} • 9> a + /i.e. , (B2) 
r.AA = { (^ 7A - - A „ 7/ ,) 75 T a A 1 '} • 9 A 0a , (B3) 

where 0, N and A denote pion, nucleon and delta fields, f n is the pion decay constant, g^, 
g-irNA and g-xAA are coupling constants, and r, M and T are matrices that couple nucleons 
and deltas into isospin 1 states, with r a r a = 3, M\M a = 2, M a M\ = 1, and T^T a = 15/4. 
For the coupling constants we use = 1-25 and the SU(A) results g^NA — ^9a/^V2 
and g nAA = 3 g A /5. We also use f n = 93 MeV, /i = 139.57 MeV, m = 938.27 MeV and 
M = 1232 MeV. 
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nucleon: 



Using the loop integrals S defined in appendix A, we obtain the following contributions 
to the nucleon scalar form factor 



a NA (r) = 2 

,2 



MA 



MttNA 



(m+M) 
QM 2 



(m+M) (2M-m) + 2^ 2 



m/i 2 (l-V 2 /2) 
(m+M) 



(B4) 



o(00) 



2m n 



[m+M) 
1 



c(oo) 



2m fi 
1 



[(m 2 -M 2 ) (m+M) (2M-m) + 2/i 2 (m 2 -M 2 ) + 6M 2 /i 2 (l-V 2 /2)] 5 



(000) 



^ 2m 
delta: 



6M 2 u 2 

(m+M) (4mM-M 2 -m 2 ) - 2^ 2 (2M-m) + j—^m i 1 -^ 2 / 2 ) 



S 



Atttt ^ 



(100) 



(B5) 



In the evaluation of the triangle diagram, the external deltas are on shell and one has the 
constraints p- u s (p) = 7- u s {p) = u s '(p')-p' = u s '{p')-^j = 0. The T matrix can be cast in 
the form 

d 4 Q 



with 



Ujy — z 



HQ-kNA 



(2^ [{Q- q /2Y-^}[{Q + q/2f-^} 
v m + M+ Q 



Kip) 



(Q + q /2r(Q-q/2y 



—9 9 ' 
p z — m z 



15 



MnAA 



2 r 



^ 2M 



p z — M z J 3 p z — m z 



(B6) 



(B7) 



(,B8) 



with p = P + Q. Performing the integrals, comparing the results with eq.(2), and going to 
configuration space, we obtain the contributions 



^A A (r) 



HQttNA 



(m + M) =,(000) =,(100) 
2M w ™ 2M ^ 
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4nU 



2M(/i 2 



o(00) _ r<( 100 ) _ c(°°°) 1 c(100) 



(B9) 
(B10) 



APPENDIX C: CHIRAL SYMMETRY 



In this appendix we show that results (B4,B5) and (B9,B10) are fully compatible with 
standard chiral power counting by means of a covariant chiral expansion [12]. It is important 



14 



to note that these expressions contain a factor (/i 3 ), which comes from the definition of the 
configuration space function S and must not be included in the counting. With this previous 
in mind, we use the following relations among integrals, 



c(oo) 



2T72 



fi 2 V 

Am 2 



c(ioo) 

T.TTTT 



2T72 



fi 2 V 

Ami 



£(000) 



2m, 



^fi-v 2 /2) + M^! 

2/im e 

r-(OO) 



2m e 

l-V 2 /2) 



c(000) 



(CI) 



2/im e 



m e m x) rn e +m x fi 



m; 



4/i 2 rrtg 



2mi 



Ami Ami 



c(ooo) 



(C2) 



sffi? = -i(i-v a /4)5£? 



Ami j 
(m 2 — m 2 ) 2 m 2 e + m 



+ 



+ 



m; 



g(100) 



(C3) 



An 2 m\ 1m 2 e Ami Am\ 

which are obtained by multiplying eqs.(A8-AlO) by P M , neglecting short range terms with a 
single pion propagator, and going to configuration space. In the case m e = M and m x = m, 
the expansion of eqs.(A8-A9) yield the leading order relations 



(ooo) (2m e /j) 

ii A7r7r 



tt(°°) 

(m 2 -M 2 ) ™ ' 



ii A7r7r — 



(2m e /i) 2 



n(oo) _ 



(m 2 -M 2 ) 2 

Truncating the expansions at 0(q A ), we find 



(C4) 
(C5) 



M/i 3 ) [(i-vV^C-^v 2 ^ 



^9ttNA 



1 2 



(M-m 



;i-v 2 /2) - 



m 

3M 2 



l-V 2 /4) 



(C6) 
S$p, (C7) 



+ 1 



(M+m) 3 



(M+m) 2 2 ^(ioo) 
4/zM 2 tM m ^™ 



M-m 



6M 



1 _ Z! ) 5 (ooo) 



Ntvtt 



^A A (r) = — 



^5Vaa 



/i(/i s 



V2 
c (100) 



2M A ™ 



_ 4 ^(000) Jj^_ r-(100) 

g A7T7T 3J\/f '^ 7r7r 



(Q 



(C9) 
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These results, except for eq. C8, which contains imaginary terms, are compatible with 
chiral counting rules. Contributions begin at 0(q 3 ) for diagrams in which internal and 
external baryons are identical and at 0(q 4 ) when this does not happen. 



APPENDIX D: CHIRAL LIMIT 



In this section we show that our model for the nucleon cx-term is consistent with the 
standard ChPT expansion. In the paper by Becher and Leutwyler[8], one finds, using our 
notation 

9g 2 A fi 3 3^ 



-4ci /i z - ~* A r " - T, (g 2 A -8c l m+c 2 m+4:C 3 m) In — 
3u 4 

aA 2 , 2 (3^_8 Cl m+4c 3 m) + 2e x , 



(Dl) 



where q and e\ are, respectively, low enegy constants (LECs) from the C/§ and la- 
grangians. The bar over e\ indicates that it has been renormalized. 
In order to expand our cr^v, we use in eq.(C4) the result 

„-2x 



c(ooo) _ 

°Ntttt — 



+ 



2x 2 mux 2 



[xK (2x) +K l (2x)} 



:i>2) 



which holds[15] for fi/m << 1. This allows integrations in eq.(10) to be performed analyti- 
cally and one finds 
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An expansion for small values of fi yields 
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where 7 is the Euler constant. This result reproduces the first three terms of eq.(Dl), 
provided one absorbs the factors proportional to \j? into the definition of ci, uses the delta 



1(3 



contributions to the Cj, which are given by 

a _ n A _ 4 9l N A m2 a _ ~ 4 91na m z\ 

Cl ~ U ' ° 2 ~ 9M2(M-m) ' ° 3 " 9(M-m) ' 1 ] 

and chooses the value R — 1/m for the cutting radius. As the renormalized constant e\ 
contains factors proportional to /i 4 , terms of this kind need not to coincide. 



ACKNOWLEDGMENTS 

The works by I.P.C. and D.O.S. were supported by CNPq and the work by G.R.S.Z., by 
FAPESP (Brazilian agencies). 



[1] T. R. Hemmert, B. R. Holstein and J. Kambor, J. Phys. G 24, 1831 (1998). 

[2] N. Kaiser, S. Gerstendorfer and W. Weise, Nucl. Phys. A 637, 395 (1998); N. Fettes and 

G-Ulf. Meissner, Nucl. Phys. A 679, 629 (2001). 
[3] V. Bernard, T. R. Hemmert and Ulf-G. Meissner, preprint hep-lat/0503022. 
[4] V. E. Lyubovitskij, Th. Gutsche, A. Faessler and E. G. Drukarev, Phys. Rev. D 63, 054026 

(2001). 

[5] H. Hellmann, Einfiihrung in die Quantenchemie (Deuticke Verlag, Leipzig, 1937); R.P. Feyn- 
man, Phys. Rev. 56, 340 (1939); S. T. Epstein, Am. J. Phys. 22, 613 (1954); J. Gasser, M. 
E. Sainio and A. Svarc, Nucl. Phys. B307, 779 (1988). 

[6] J. Gasser, H. Leutwyler and ME. Sainio, Phys. Lett. B 253, 252 (1991); 253, 260 (1991). 

[7] J. Gasser, M.E. Sainio and A. Svarc, Nucl. Phys. B 307, 779 (1988). 

[8] T. Becher and H. Leutwyler, Eur. Phys. Journal C 9, 643 (1999); JHEP 106, 17 (2001). 

[9] M.R. Robilotta, Phys. Rev. C 63, 044004 (2001). 
[10] A.I. L'vov, S. Scherer, B. Pasquini, C. Unkmeir and D. Drechsel, Phys. Rev. C 64, 015203 
(2001). 

[11] T. H. R. Skyrme, Proc. Roy. Soc. London A247, 260 (1958); A252, 236 (1959); A260, 127 

(1961); Nucl. Phys. 31, 556 (1962). 
[12] R. Higa and M. R. Robilotta, Phys. Rev. C 68, 024004 (2003). 
[13] C. Fronsdal, N. Cim. Suppl. 9, 416 (1958). 



17 



[14] G. Hohler, group I, vol.9, subvol.b, part 2 of Landolt-Bornstein Numerical data and Functional 

Relationships in Science and Technology, ed. H.Schopper (1983). 
[15] R. Higa, M. R. Robilotta and C. A da Rocha, Phys. Rev. C 69, 034009 (2004). 



18 



